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Abstract — Stochastic network calculus requires special care 
in the search of proper stochastic traffic arrival models and 
stochastic service models. Tradeoff must be considered between 
the feasibility for the analysis of performance bounds, the use- 
fulness of performance bounds, and the ease of their numerical 
calculation. In theory, transform between different traffic arrival 
models and transform between different service models are 
possible. Nevertheless, the impact of the model transform on 
performance bounds has not been thoroughly investigated. This 
paper is to investigate the effect of the model transform and to 
provide practical guidance in the model selection in stochastic 
network calculus. 

Index Terms — Stochastic Network Calculus, Model Transform, 
Performance 

I. Introduction 

Performance has always been one of the major concerns 
in networking systems. Mathematical models for quantitative 
evaluation of network performance, however, have remained 
as a slow-paced research area. A.K. Erlang published the first 
paper on queuing theory in 1909 [1], and since then queuing 
theory has been developed and applied in a wide variety 
of applications. In particular, it has been the foundation in 
performance modeling and evaluation of telecommunication 
systems and has been applied broadly in the performance 
analysis of computer networks. Nevertheless, with the advance 
of the Internet technology, the assumptions behind the tractable 
queuing models may not hold anymore. Despite the research 
efforts in the last one hundred years, the tractable models with 
traditional queueing theory consist of only a minority of practi- 
cal network problems. The research community is in dire need 
of new mathematical models for network-wide performance 
evaluation where the Markovian property in traffic arrivals or 
services may not hold. 

Network calculus is one of such new analytical techniques. 
The theory of network calculus was introduced in early 
1990s [5 1 for network performance evaluation. Unlike the 
traditional queueing theory which aims at obtaining exact 
analytical results, network calculus focuses on the analysis 
of performance bounds using the cumulative amount of traffic 
arrivals or services. Since network calculus usually does not 
assume particular distributions on traffic arrivals or service 
times, it can obtain broadly applicable performance results. 
Network calculus has evolved along two tracks- determinis- 
tic El, E6] and stochastic flU, 03, 03, EQl- The determin- 
istic network calculus is to obtain the worst case performance 
bounds, which may be too loose to be useful in practice. Due 
to this reason, research on this direction gradually fades out. 
To overcome the problem, stochastic network calculus was 
developed. Nevertheless, due to some special difficulties fTOll . 
03, basic properties of stochastic network calculus have been 



proved only in recent years flU, iflOll . ifTTIl . Although the major 
theoretical barriers have been cleared, it is still unclear whether 
or not stochastic network calculus will be broadly accepted by 
network practitioners. 

Without the driving force from real applications, broad ac- 
ceptance of stochastic network calculus as a valuable technique 
for performance evaluation may not be optimistic. One of the 
major practical challenges is the lack of effective algorithms 
to calculate and compare the performance bounds. After all, 
what really matter to network engineers are the meaningful 
numerical results instead of the complex equations. Although 
there are some efforts using Legendre transform Q, @ 
to simplify the calculation of major operations in network 
calculus, the treatments are far from sufficient to tackle the 
difficulties in the stochastic network calculus, where we are 
often faced with multiple tradeoffs. 

Specifically, three tradeoffs must be considered in stochastic 
network calculus. First, there is a tradeoff between the sim- 
plicity of deriving performance bounds and the difficulty in the 
numerical calculation of the bounds. It is known that in order 
to derive performance bounds easily, we need to put extra con- 
straints on the traffic and the service models [ 1 1 1. For instance, 
we may need to put some constraints in the cumulative traffic 
arrivals/services, e.g, we may change the calculation from the 
form of Prob{ft > 0} to the form of Prob{sup t f t > 0}, 
which is usually not equal to sup t Prob{f t > 0}. Note 
that sup is the supremum (i.e., least upper bound) operation. 
Prob{sup t ft > 0} thus represents an instantaneous property 
and is generally hard to calculate. Second, there is a tradeoff 
between the usefulness of the traffic (service) models and the 
hardness of searching for these models in real applications. 
This tradeoff is closely related to the first one. In general, 
it is easy to obtain the traffic model (or service model) 
directly from the distribution of packet inter-arrival times 
(or the distribution of the service times). Introducing extra 
constraints on the traffic arrival or service model, e.g., adding 
the sup operation in the model iflOl . requires that we either 
perform model transform ifTTI or search for the models using 
queueing analysis methods [U, 03- Third, we must consider 
the tightness of performance bounds in a stochastic sense. In 
stochastic network calculus, we need to weigh a performance 
bound regarding its tightness and its bounding function, e.g., 
we need to avoid poor claims like "the probability that the 
delay is larger than 30 seconds is less than 90%," which is 
not helpful in practice. 

Handling the above problems has been a very tricky and 
intimidating task. Without a clear guideline, it may not be easy 
to use stochastic network calculus in real-world problems. We 
are thus motivated in this paper to analyze the above tradeoffs 
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and provide clear guidance on the tricky model selection 
and model transforms. Although there is a software package, 
called DISCO fl9l , to "automatically" derive deterministic 
performance bounds once model parameters are given, such a 
software tool has not been seen so far for stochastic network 
calculus, due to the above tricky tradeoffs. We expect this 
paper can also clear the road for people who intend to build 
a software package for stochastic network calculus. 

II. Related Work 

As a new theory for performance evaluation, network cal- 
culus has been developed along two tracks: deterministic and 
stochastic. Deterministic network calculus |2j, lfl6ll is to search 
for the worst-case performance bounds, which in many cases 
are too loose to be useful. Stochastic network calculus Q, 
1 10 1, [11] tries to derive tighter performance bounds, but with 
a small probability the bounds may not hold true. In practice, 
the bounds obtained via stochastic network calculus may be 
more useful, since such bounds present network engineers with 
a mechanism to utilize statistical multiplexing gain. 

It seems a strange phenomenon that most papers on stochas- 
tic network calculus mainly focus on the theoretical devel- 
opment. Various types of calculus are proposed to analyze 
the performance bounds in a stochastic sense |4), Q, |6|, 
ifTTl . 1171 . Different approaches have been used, for example, 
the effective bandwidth H, flTl . moment generating func- 
tions [6], Martingale inequality [12]. In [10], a stochastic 
network calculus is built with quite generic abstract traffic 
models and service models. Each calculus, without an ex- 
ception, has been demonstrated to be effective and useful 
for some given application scenarios. The quest on using 
stochastic network calculus to solve queueing problems has 
been remaining active Q, fl2l . In contrast, the applications 
of stochastic network calculus were left behind. It is abnormal 
that the papers on theoretical development outnumber the ones 
on realistic applications of this theory. 

The call for a guidance on building suitable stochastic traffic 
and service models that are simple to obtain and easy to 
calculate remains unanswered. While many research efforts 
are being devoted to obtaining tight stochastic bounds close to 
the exact solutions for special cases [3], we in this paper divert 
to pursuing the simplicity of the model building methodology. 

III. Background of Stochastic Network Calculus 

A. Notation 

We first introduce the notation and key concepts of stochas- 
tic network calculus 1101 . IfTTl . IfTTl . Throughout this paper, 
we assume that all arrival curves and service curves are non- 
negative and wide-sense increasing functions. Conventionally, 
A(t) and A*(t) are used to denote the cumulative traffic that 
arrives and departures in time interval (0,i], respectively, and 
S(t) is used to denote the cumulative amount of service pro- 
vided by the system in time interval (0, t]. For any < s < t, 
let A(s,t) ee A(t) - A(s),A*(s,t) = A*(t) - A*(s), and 
S(s,t) ee S(t)-S(s). By default, A(0) = A*(0) = 5(0) = 0. 

We denote by T the set of non-negative wide-sense increas- 
ing functions, i.e., 



and by T the set of non-negative wide-sense decreasing 
functions, i.e., 

F = {f(-):W<x<y,0<f(y)<f(x)}. 

For any random variable X, its distribution function, de- 
noted by 

F x (x) ee Prob{X < x}, 

belongs to T, and its complementary distribution function, 
denoted by 

F x {x) ee Prob{X > x}, 

belongs to T. 

During model transform, we may put a stronger requirement 
on the bounding function. We denote by Q the set of functions 
in T where for each function g(-) € Q, its nth-fold integration 
is bounded for any x > and still belongs to Q for any n > 0, 
i.e., 

/ /-oo \ n 
3 = {<?(•) : Vn > 0, ( J dyj g{y)eQ}. 

B. Operators 

The following operations defined under the (min, +) alge- 
bra Q, 0, lfl6l will be used in this paper: 

• The (min, +) convolution of functions / and g is 

{f®9){t)=mi{f{s)+g{t-s)}. (1) 

0<s<t 

• The (min, +) deconvolution of functions / and g is 

(/0s)(*)=sup{/(* + «)-<?(«)}. (2) 

In addition, we adopt: 

• [x] + = max{x, 0}, 
« [x]i = min{x, 1}. 

C. Performance Measures, Traffic and Server Models 

The following measures are of interest in service guarantee 
analysis under network calculus: 

« The backlog B(t) in the system at time t is defined as: 

B(t) = A(t)-A*(t). (3) 

• The delay D(t) at time t is defined as: 

D(t) = inf{r > : A(t) < A*(t +r)}. (4) 

Stochastic traffic arrival curve and stochastic service curve 
are core concepts in stochastic network calculus with the for- 
mer for traffic modeling and the latter for server modeling. It is 
worth noting that the deterministic traffic arrival curve and the 
deterministic service curve under the (deterministic) network 
calculus are a special case of their corresponding stochastic 
definition. In the literature, there are different definitions of 
stochastic arrival curve and stochastic service curve flOl . ifTTIl . 
For traffic arrival models, we have: 

Definition 1: The t.a.c. model: A flow A(t) is said to 
have a tjaffic-amount-centric (t.a.c.) stochastic arrival curve 
a G J- with bounding function / 6 T, denoted by 



^ = {f(-) < x < y,0 < f(x) < f(y)}, 
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if for all t > s > and all x > 0, it holds Q0), 02 

Pro6{yl(s, t) - a(t - s) > x} < f(x). (5) 

Definition 2: The v.b.c. model: A flow A(t) is said to 
have a virtual-backlog-centric (v.b.c.) stochastic arrival curve 
a £ J 7 with bounding function / G T, denoted by 

A ~„6< f,a >, 
if for all t > and all x > 0, it holds OH, fTTl 

Prob{ sup [-4(s,t) - a(f - s)] > x} < /(x). (6) 

0<s<t 

Definition 3: The m.b.c. model: A flow A(t) is said 
to have a maximum-virtual-backlog-centric (m.b.c.) stochastic 
arrival curve a £ J 7 with bounding function / G T, denoted 
by 

A ~ m j< /, a >, 
if for all t > and all x > 0, it holds OH), 02 

P?'ofo{ sup sup [A(u, s) — a(s — u)] > x} < f(x). (7) 

0<S<(0<H<S 

For service models, we have the followings. 

Definition 4: The w.s. model: A server is said to provide 
a flow A(t) with a weak stochastic (w.s.) service curve /3 G T 
with bounding function g G .F, denoted by 

«5 ~ u, s < 3, /3 >, 

if for all i > and all x > 0, it holds OH, 02 

Prob{A ® /3(f) - A*(t)] >x}< g(x). (8) 

Definition 5: The s.c. model: A server is said to provide 
a flow A(t) with a stochastic service curve (s.c.) (3 G T with 
bounding function g G .F, denoted by 

S ~ sc < 9,(3 >, 
if for all f > and all x > 0, it holds 021, 02 

Prob{ sup [A <gi P(s) - A*(s)] > x} < ff (x). (9) 

0<s<t 

Definition 6: The s.s.c. model: A server is said to to 
provide a strict stochastic service curve (s.s.c.) (3 E J- with 
bounding function g G F, denoted by 
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if during any period (s,t] the amount of service S(s,t) 
provided by the server satisfies ifTUl . iTTTl 



Prob{S(s,t) < P(t -s)-x}< g(x). 



(10) 



With the above definitions, various properties of stochastic 
network calculus, including the stochastic backlog bound and 
the stochastic delay bound, have been proved (e.g., see |10|, 

02,03). 

Several natural questions arise: Why should we need dif- 
ferent forms of traffic arrival models and service models? 
Can a traffic (service) model be transformed to another traffic 
(service) model? What is the impact of model transform on 
performance analysis? 



The first two questions have been answered in [10], ifTTIl . 
Briefly speaking, some models are too weak to be useful 
in the performance bound analysis. For instance, from the 
t.a.c. model it is hard to obtain the stochastic backlog 
bound, because according to Lindley equation [14|, B(t) = 
sup 0<s<t {A(s, t) — S(s,t)}, which requires the calculation 
of sup. The value of sup is not readily obtainable from the 
t.a.c. model. As such, we may put more constraints on the 
traffic model, such as those in the v.b.c. model and the m.b.c. 
model. Regarding the second question, it has been shown that 
different models can be transformed to each other with the 
theorems introduced in 1111 . The last question, however, has 
not been touched and is the main focus of the rest of the paper. 

IV. Tradeoffs in Model Transform 

In this section, we illustrate the tradeoffs in the selection 
of a proper model. We start from the t.a.c traffic arrival 
model and the s.s.c. service model, because both of them 
have the most intuitive meaning and can be obtained easily 
from the distribution of packet inter-arrival times of the input 
flow and the distribution of packet service times of the server, 
respectively. 

We shall ignore the transform from m.b.c — > v.b.c. t.a.c. 
and the transform from s.c — > w.s., because a stronger 
model (i.e., a model with more constraints) implies a weaker 
model ifTTl . 

A. The Transform from t.a.c. to v.b.c. 

Although the t.a.c. model is the most simple model, it is 
not easy to derive performance bounds with this model IfTTl . 
As such, we need to transform it to a stronger model, e.g., the 
v.b.c. model or the m.b.c. model. 

Remark 1: It is very likely that after the transform the 
values of the bounding function go to a large probability 
value or even 1 when time t goes to oo, especially when the 
bounding function of the t.a.c. curve is dependent on time t. 
This is the so-called time-increasing problem on the bounding 
function. To avoid this problem, it is suggested [ 17 1 that there 
should exist a time scale T enforced on the traffic and the 
service. We follow same idea assume that all traffic arrival 
curves and service curves are enforced on the maximum time 
scale T, e.g., the constraint of v.b.c. traffic curve becomes 

Prob{ sup [A(s, t) - a(t - s)] > x} < f(x). 

t-T<s<t 

Later, we will discuss other methods to tackle this problem. 

Lemma 1: If a stationary traffic flow has a t.a.c stochastic 
arrival curve a G J- with bounding function / G Q, it also 
has a v.b.c. stochastic arrival curve ae G T with bounding 
function f G Q, where for any 9 > 



ag = a(t) + 6 ■ t, 



x+T8 



f(y)dy 



(ID 

(12) 
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Proof: Since the traffic flow is stationary, for any 9 > 
0, ti >T,t< T, we have 

sup {^4(s, £) — ag(t — s)} 

t x -t<s<ti 

= sup {A(s, t) — a$(t — s)} 

0<s<t 

< sup {A(s, t) — Otg(t — s)} 

0<s<T 

< sup {A(s,t) -a e (t-s)} + . (13) 

0<s<T 

Since for any x > 0, Prob{[A(s,t) — a e (t — s)] + > x} = 
Prob{A(s, t) - a e (t - s) > x} < f(x + 6 • (t - s)), we have 

Pro6{ sup i) -ag(t- s)} > x} 

ti-t<s<ti 

<Prob{ sup {A(s,T) - a e {T - s)} > x} 

0<s<T 

T 

<J2 Prob i[ A ( s , T ) ~ a e {T - s)] + > x} 

T 

<J2f(x + e-(T- s )) 



s=0 



< 



x+TB 



f{y)dy. 



(14) 



The theorem holds since the probability has to be smaller than 
1. ■ 

Essentially, Lemma [TJ indicates that for a flow following a 
t.a.c. traffic arrival curve, we can model the same flow with 
a series of v.b.c. traffic arrival curves. An interesting question 
is that among these v.b.c. traffic arrival curves, which one is 
the best with which we can obtain the tightest performance 
bounds? To evaluate, we need to formally define the tightness 
of traffic arrival curves in the stochastic sense. 

Definition 7: (Stochastic tightness of traffic arrival 
curves) Assume that a traffic flow A(t) follows a t.a.c. (or 
v.b.c, m.b.c.) traffic arrival curve ol\ with bounding function 
/i as well as a t.a.c. (or v.b.c, m.b.c, respectively) traffic 
arrival curve a-i with bounding function f%. We call the curve 
oi\ is stochastically tighter than the curve ct2 within a tolerance 
bound e > 0, denoted by ct\ < e a.2, if for all t > and all 
x > 0, there hold 



e*i(t) 



<a 2 (t) if i = 
< a2(t) otherwise, 



and 



A(x) <h(x) + e. 



(15) 



(16) 



If e = 0, we also say that ct\ is absolutely tighter than a-i- 

Generally speaking, we need to make a tradeoff between the 
arrival curve and its bounding function. From Lemma [Tj we 
observe that the series of v.b.c. traffic arrival curves depends 
on the value of 9. We should not select a very loose traffic 
arrival curve (i.e., a very large 9 value) to make the bounding 
function small; nor should we use a very tight traffic arrival 
curve (i.e., a very small 9 value) such that the bounding 
function becomes not useful. For example, "the probability 
that a certain event occurs is no larger than 1" is meaningless. 



We have the following theorem to determine another tighter 
traffic arrival curve based on an existing traffic arrival curve 
and a given acceptable range on the bounding function. 

Lemma 2: Assume that a traffic flow has a t.a.c stochastic 
arrival curve a G T with the bounding function / e Q. 
Assume that one of its corresponding v.b.c. stochastic arrival 
curves is ag 1 £ J- with the bounding function f 01 £ Q. We can 
model it with another v.b.c. stochastic arrival curve ag 2 £ T 
with the bounding function f 8 ' 2 £ Q such that ag 2 < e ae 1 if 
there exists 62 < 9\ for any x > satisfying 

X+T82 ^ i-x+TBx ^ 

T f(y)dy - / T f(y)dy < e. (17) 

: 02 Jx "1 

Lemma [2] is easy to prove based on Lemma Q] and Defini- 
tion [7] 

Remark 2: Based on Lemma if no such 6*2 could be 
found, we call the arrival curve ag 1 the stochastically tightest 
within the tolerance bound e. Given any v.b.c. traffic curve, 
Lemma |2] is useful in searching for a tighter stochastic arrival 
curve, if exists. For example, assume that we have a v.b.c 
traffic curve < f 01 , a$ 1 >, where d\ is known. Denote 



2, 



X+T82 



-K-f{y)dy 

t>2 



x+TOt 



g^f(y)dy- 



Setting a lower threshold value on x, say x, and a tolerance 
bound e, we can check if the equation </>(#2, n) — e nas a pos- 
itive root on 62- If no solution could be found, < f 01 ,ctg 1 > 
is the tightest v.b.c. curve within the tolerance bound e. 

B. The Transform from v.b.c. to m.b.c. 

Lemma 3: If a traffic flow has a v.b.c stochastic arrival 
curve a £ T with bounding function / £ Q, it also has a 
m.b.c. stochastic arrival curve ag with bounding function / £ 
Q, where for any 9 > 



ag = a(t) + 9 ■ t, 



fix) 



f(y)dy 



-8T 



(18) 
(19) 



The proof of Lemma [3] is similar to that of Lemma [TJ 
and is omitted. Due to the similarity between Lemma [3] and 
Lemma [TJ we can slightly revise Lemma [2] so that we can 
check whether there exists a stochastically tighter m.b.c. curve, 
given an existing m.b.c. curve. 

C. The Transforms from s.s.c. to w.s and s.c 

Lemma 4: Consider a server that provides a stochastic 
strict service curve f3 £ T with bounding function g(x) £ T . 

1) The server also provides a weak stochastic service curve 
j3(t) with the same bounding function g(x). 

2) If g{x) £ Q and the input and output processes are both 
stationary, the server provides a stochastic service curve 
j3-g with bounding function g {x), where for any 9 > 0, 



/3_ e = p{t) -9-t, 



f(y)dy 



-0-T-\ 



(20) 
(21) 



Proof: Please refer to chapter 4 of ifTTIl for the proof of 
the first part. 
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For the second part, we first have for any t > s, 

A®/3(s) <A®l3(s)-6(t-s), 
and hence for any t\ > T,t < T, 

Prob{ sup [A ® fa-e(s) - A*(s)} > x} 

t 1 -T<s<t 1 

<Prob{ sup [A <8> /3_ e (s) - A*(s)} > x} 

0<s<T 

<Prob{ sup [A /3(s) - A*(s) - 9 • s]+ > a; - • T} 

l<s<T 

T 

<^Pro6{[A®/3(s)-A*(s)-0-s]+ > x - 9 ■ T} 



<J2f( x - e - T + 9 - s ) 

s=l 

1 r 

/ f(y)dy 

Jx-e-T+6 



(22) 
(23) 



Since the probability cannot be larger than 1, the second part 
is proved. Note that the inequality d22b is true due to the first 
part of the theorem. ■ 
Lemma [4] indicates that if we transform the s.s.c service 
model to the s.c. service model, we obtain a series of s.c. 
curves and are faced with the problem of selecting a "good" 
s.c. curve for performance analysis. Similar to the transform 
of traffic models, we need to define the tightness of service 
curves. 

Definition 8: (Stochastic tightness of service curves) 

Assume that a service provided by a system follows an s.s.c. 
(or w.s., s.c.) service curve fa with bounding function g\ 
as well as an s.s.c. (or w.s., s.c, correspondingly) service 
curve fa with bounding function g 2 . We call the curve fa 
is stochastically tighter than the curve fa within a tolerance 
bound e > 0, denoted by fa > c fa, if for all t > and all 
x > 0, there hold 



>fa{t) if t = 
> fait) otherwise, 



and 



9i(x) < 92{x) + e. 



(24) 



(25) 



If e = 0, we also call that fa is absolutely tighter than fa. 

When we select a good stochastic service curve, we again 
need to make the tradeoff between the tightness of the service 
curve and the usefulness of the bounding function. Similar 
to the previous section, we have the following theorem to 
help select a good stochastic service curve after the model 
transform. 

Lemma 5: Assume that a server provides an s.s.c service 
curve (3 E T with the bounding function g E Q. Assume that 
one of its corresponding s.c. service curves is /3-$ 1 E T with 
the bounding function g 01 E T. We can model it with another 
s.c. service curve fa-o 2 E J- with the bounding function g 02 E 
F such that f3-g 2 > e fa-g 1 if there exists 6 2 < 9\ for any 
x > satisfying 

f x 1 , . F 1 



f(y)dy 



ix-9 2 +e 2 



-f(y)dy < e. (26) 



X— 01+01 



Lemma [5] is easy to prove based on Lemma [4] and Defini- 
tion m 

V. The Impact of Model Transform on 
Performance Evaluation 

In this section, we illustrate the impact of using different 
traffic arrival/ service curves on the performance evaluation. 
We only use the output characteristics and the service guar- 
antee (e.g., stochastic bounds on delay and backlog) as the 
examples. The impact on other properties such as the leftover 
services is omitted to save space. 

Lemma 6: 1) If for any x > 0, fi(x) < fzfa) + e\ and 
9i(x) < g2(x) + e 2 , then /1 <g> g x {x) < fx ® gi(x) + 
ei + £2- 

2) If for any x > 0, ctx{x) < a 2 (x) and fa(x) > fa{x), 
then a\ fa{x) < a 2 fa{x). 

3) If for any x > 0, a\{x) < a%{x) and fa(x) > fa{x), 
then h(ai,fa) < h(a 2 ,fa), where h(a, fa) represents 
the maximum horizontal distance between functions a 
and fa i.e., 

h(a, fa = sup{inf{r > : a(s) < fas + t}}. 

Proof: For the first part, we have: 

h ® 9i(x) = w£ {fi{y)+gi(x-y)} 

0<y<x 

- ni n t y^y) + e i + 92& -y) + £2} 

0<y<x 

= n3 n t {h(y) + 92{x - y)} + €1 + e 2 

0<y<x 
= I2 ®92{x) +61 +£2- 
For the second part, we have: 

a x fa{x) = sup{ai(x +y)- fa{y)} 
y>o 

< sup{a 2 (x + y) - fa(y)} 
y>o 

= a 2 fa(x). 
To prove the last part, we have: 

h{ai,fa) = sup{inf{r > : ai(s) < fa(s + t}} 

s>0 

< sup{inf{r > : a 2 (s) < fa(s + t}} 

s>0 

< sup{inf{r > : a 2 (s) < fa(s + t}} 

s>0 

= h(a 2 ,fa). 



From Lemma|6]and Theorem 5.12 in ifTTI . it is easy to have 
the following theorem: 

Lemma 7: Consider a system with input A. Assume that 
A ~t>6< f\,ot\ > as well as A ^ v b< f2,&2 >, and the 
system provides the input A with a service that can be modeled 
by S ~ sc < g 1 ,fa > and S - sc < g 2 ,fa >. If a x is 
stochastically tighter than a 2 within the tolerance bound t\ 
and fa is stochastically tighter than fa within the tolerance 
bound € 2 , then the output A* can be modeled by: 

1) A* ~„ 6 < /1 g u ax fa >, or 

2) A* ~ vb< f l (g, g 2 , ai fa >, or 
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3) A* ~ vb < / 2 <8> #1, a 2 0i >, or 

4) A* ~ vb < f 2 o 92i a x >, 

among which A* ~ vb < /i0gi, ai0/3i > is the stochastically 
tightest within tolerance bound e\ + e 2 . 

From Lemma [6] and Theorem 5.4 in IfTTI . it is easy to have: 
Lemma 8: Consider a system with input A. Assume that 
A r ^ v b< fiiO-i > as well as A ^ vb < f2i a 2 >> an d the 
system provides the input A with a service that can be modeled 
by S ~ ws < g l ,0i > and S ~ ws < g 2 ,/3 2 >■ If «i is 
stochastically tighter than a 2 within the tolerance bound ei 
and (3i is stochastically tighter than /3 2 within the tolerance 
bound e 2 , then for all £ > and x > the delay D(t) can be 
bounded by : 

1) Prob{D(t) > h( ai + x, ft)} < / 1 ^(x), or 

2) Prob{D(t) > h{a x + x, f3 2 )} < f 1 5 2 (x), or 

3) Prob{D(t) > h(a 2 + x, 0i)} < f 2 g^x), or 

4) Profe{D(i) > h(a 2 + x, #,)} < f 2 5 2 (a;), 

among which Prob{D(t) > ft,(ai + < / 1 g 1 ^) is 

the stochastically tightest bound in the sense that h(ai+x, (3i) 
is the smallest for all h(ai + x, f3j),i, j — 1,2, and f 1 g 1 — 
e i — e 2 is also the smallest for all /* g\ i,j = 1, 2. 

Regarding the backlog bound, we have: 

Lemma 9: Consider a system with input A. Assume that 
A ~ v b< fi,o.i > as well as A ^ vb < / 2 ,a 2 >, and the 
system provides the input A with a service that can be modeled 
by S ~ ws < g x ,0i > and S ~ ws < .g 2 ,/3 2 >■ If «i is 
stochastically tighter than a 2 within the tolerance bound e\ 
and (3i is stochastically tighter than /3 2 within the tolerance 
bound e 2 , then for all t > and x > the backlog B(t) can 
be bounded by : 

1) Prob{B(t) >x}<f 1 ®g 1 {x-a 1 (Z) /?i(0)), or 

2) Prob{B(t) > x} < f 1 g 2 \x - a x /3 2 (0)), or 

3) Prob{B(t) >x}<f®g 1 {x-a 2 <3 /?i(0)), or 

4) Prob{B(t) >x}<f 2 ® g 2 (x - a 2 /3 2 (0)), 
among which Prob{B(t) >i}</ 1 ®ff 1 (i-ai0 /MO)) is 
the stochastically tightest bound in the sense that f 1 ®g 1 (x — 
ai /3i(0)) is the smallest for all p g' J (x -a{Q) ft (0)) - 
ei - £2,«, j = 1,2. 

Lemma [9] can be easily proved based on Theorem 5.1 
in IfTTI . Lemma|6] and the fact that if /, g 6 T then f®g G T . 

VI. Integrating Queueing Theory and Stochastic 
Network Calculus 

In j3|, [12], performance bounds close to the exact so- 
lutions with traditional queueing theory have been derived 
with stochastic network calculus. These methods derive perfor- 
mance bounds based on concentration inequality and martin- 
gale inequality. In addition, it is required to solve optimization 
problems in the selection of the best parameters to obtain tight 
bounds 0. These methods are effective but the complexity in 
model buildup may make them not easily accessible to a broad 
range of users. 

We suggest that model building and performance analysis 
should be decoupled to ease the application of stochastic 
network calculus. More specifically, it should be a two-phase 
process. In the first phase, tight traffic arrival and service 
models should be built, and all the complexity requiring the 
help of traditional queueing analysis should be dealt with in 



this phase. The second phase is to automatically calculate 
performance bounds with standard operations in stochastic 
network calculus. It should be feasible to build a software 
tool to handle the bound analysis in the second phase. 

Traditional queueing theory plays an important role in the 
first phase. While it is hard to suggest a general approach 
suitable for all applications, the following guidelines should 
be helpful. 

• Rule 5: Regarding traffic arrival curves, we can build a 
virtual queueing system, with the traffic arrival process as 
the input and the arrival curve as the service of the virtual 
queueing system. Building the v.b.c (or m.b.c.) model 
is equivalent to finding the queue length (or maximum 
queue length, respectively) of the virtual queueing system 
with respect to time t. This is a traditional queuing theory 
problem, and tight bound can usually be found. 

• Rule 6: Regarding service curves, both w.s. and s.c. 
models are coupled with the input process (in definition). 
As such, we suggest starting with the s.s.c. model, which 
has a straightforward meaning. From the s.s.c. model, 
we can directly get the w.s. model. In the following, we 
propose another new service model, which is easy to find 
with queuing theory and can be used to obtain both w.s. 
and s.c. models directly. 

Definition 9: (Virtual backlog stochastic strict service 
curve). A system is said to be a virtual backlog stochastic strict 
server providing service curve /3(t) with bounding function 
g(x) G T, denoted by S ~ v bssc< g(x),f3(t) >, if during 
any period (s, t] the amount of service S(s, t) provided by the 
system satisfies 

Prob{ sup \p{t - s) - S{s, £)] > x} < g(x), (27) 

0<s<t 

for any x > 0. 

We have the following theorem: 

Theorem 1: Consider a system that is a virtual backlog 
stochastic strict server providing service curve f3(t) with 
bounding function g(x) G T. It also provides an s.c. stochastic 
service curve j3(t) with the same bounding function g(x). 
Proof: For any t < 0, V0 < s < t, we have 

• Case 1: s is not within any backlogged period. In this 
case, A*(s) — A(s). Therefore, 

A p(s) - A*(s) < A(s) + f3(0) - A*(s) < 0. (28) 

The last inequality is because (3(0) needs to lower bound 
the service, which is S(0) = by default. 

• Case 2: s is within a backlogged period. Assume that the 
backlogged period starts from so < s. Then A*(sq) = 
A(so), and 

A0/3(s) - A*(s) 
<A(s )+(3{s~s )-A*(s) 
= p(s-s ) + A*(s )-A*(s) 
= 0(s-so)-S(s o ,s) (29) 

Combining both cases, we have 

A 0(a) - A*(s) < I3{s - Sl ) - S(sx,s) 
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holds for all s G [0,t] and < si < s. Therefore, 

[/3( s - Sl )-5(s 1)S )] 



sup [A<g>/3(s) - A*(s)] 

0<s<t 



< sup 

0<s<t,0<si<s 

= sup [/3{t - s) - S(s,t)], 

0<s<t 



(30) 



from which and the definitions of s.c. and v.b.s.s.c. service 
curves the theorem is proved. ■ 
Regarding service curves, we have the following suggestion 
based on Theorem [T] 

• Rule 7: For a server providing service process S(t), we 
can build a virtual queueing system, with the service 
curve f3(t) as the input and the service process S(t) as 
the service of the virtual queueing system. Building the 
v.b.s.s.c. model is equivalent to finding the queue length 
of the virtual queueing system with respect to time t. 
Again, this is a traditional queueing theory problem. With 
the v.b.s.s.c. model, we can immediately obtain the s.c. 
or w.s. model. 

VII. An Example and Guidelines 

We use the typical M/M/l queueing model as an example 
to illustrate how an intuitive and nature model transform 
may result in poor performance bounds. Drawn from this 
example, we summarize the implications of model transform, 
and provide the insights on the possible problems in the model 
selection and transform. 

A. An Example 

Assume that a flow has fixed unit packet size. Assume that 
its packets arrive according to a Poisson arrival process with 
mean rate A. Assume that the service time of each packet in 
the server has an exponential distribution with mean i, where 
H > A so that the system is stable. 

This simple M/M/l model has exact solutions on per- 
formance bounds with traditional queueing analysis [14] and 
similar results can be obtained with special treatment on 
traffic models and service models of stochastic network cal- 
culus J3). Nevertheless, we would start from some simple and 
"intuitively" right models and are interested in the impact 
of model transform on performance bounds. Naturally, the 
Poisson traffic arrival process and the exponential server can 
be easily modeled by a t.a.c. traffic arrival curve and an s.s.c. 
service curve. To obtain the performance bounds, we then need 
to transfer the t.a.c. curve to a v.b.c. curve, according to the 
theorems in the previous section. 

Traffic Model: It is easy to see that in any time interval 
(s, s + t], for any x > 0: 



Prob{A(t) 



Xt>x}< 

k=\x+\i\ 



-At 



(Xt) k 



kl 



(31) 



It is hard to transform the above t.a.c. traffic arrival curve to 
other traffic models due to the sum in the bounding function. 
To make the model transform easy, we derive a simpler 
bounding function based on the Poisson approximation ifTSl - 



Lemma 10: A Poisson arrival process A{i) with mean rate 
A is bounded by 



Prob{A{t) -Xt> x} < e ^-(^+x)in^± 



(32) 



The proof of LemmafTOlis in the appendix. Note that when t 
goes to oo, the right-hand side of d32l goes to 1, indicating that 
we have the time-increasing problem with the above bounding 
function. To avoid the problem, we assume that there is a 
maximum time scale T enforced on the traffic arrivals and the 
service. Therefore, for any s > and any < t < T, 

Prob{A(s, s + t)~Xt>x}< e x-(XT+x)in^sl _ (33) 

With Theorem Q] we can find a series of corresponding 
v.b.c. curves A(t) ^ v b< ctg, fg >, where for any < t < T 



a e = (X + 9)- 1, 



x+8T 



e y- (XT+ y)In^pl dy]i 



(34) 
(35) 



Service Model: During any backlogged period (s, s + t], 
it is known that the packets departing from the server has an 
exponentially distributed inter-arrival times with mean —, that 
is, the departure has a Poisson process with mean - during 
the backlogged period. Therefore, 



Prob{S(s, s + t) =n} 
from which we can get for j3(t) = fit, 



\i n t n 



-fit 



(36) 



\lit-x] k k 

Prob{S(s, s + t)< f3(t) - x) < J2 ^- e_ "*- ( 37 ) 



k=0 



k\ 



Using Poisson approximation, the bounding function g(x) 
becomes: 



g(x) = 1 — e ^ ' f ,T 



(38) 



If a server has an s.s.c model, it also has the s.c. mode 
with the same service curve and bounding function. Having 
the combination of a v.b.c traffic model and an s.c. model, we 
can then derive performance bounds. 

Figs 03 and |2] show the numerical results on delay bounds 
for an M/M/l queueing system, where the average traffic 
arrival rate A is equal to 20 packets per second and the average 
service rate /i is equal to 25 packets. The results demonstrate 
that the bounding function is very sensitive to the maximum 
time scale T. When T is large, e.g, 4 seconds, the bounding 
function becomes very loose. Comparing Fig[T]and Fig. [2] we 
note that the impact of 9 on the delay bounds is not significant, 
since the value of 9 must be bounded by /i — A so that the 
delay bound does not escape to infinity. 

We can use traditional queueing analysis method to obtain 
the exact solution for the M/M/l model. From [14|, the 
packet delay in the M/M/l system, D(t), is an exponential 
random variable with mean -r^zj, i.e., Prob{D{t) > x} = 
6 -(m--M. This is the exact solution, as it gives the exact 
distribution of packet delay. From Figs [T] and [21 it is easy to 



s 





Fig. 1. The impact of T on the bounding function (9 = 0.1) Fig. 2. The impact of T on the bounding function (8 = 0.8) 



see that the exact bound obtained with the traditional queueing 
analysis method is much tighter than that with stochastic 
network calculus, if the stochastic traffic arrival and service 
models are not selected properly as above. 

B. Guidelines 

From our analysis and numerical example, we summarize 
the following practical guidelines for model selection and 
transform in stochastic network calculus. 

• Rule 1: In general, the transform from a weak model 
(e.g., the t.a.c. model ) to a strong model (e.g., the v.b.c. 
model) will result in a looser bounding function. This 
problem cannot be solved by optimizing the 9 value in the 
model transform. If a strong model can be found directly 
with some methods like those in ||8), ifTTI . |[T3l . never 
rely on model transform to obtain the strong model from 
a weak model. 

• Rule 2: When the bounding functions are time dependent, 
it is likely that the time-increasing problem will occur. 
The method of limiting the time with a maximum time 
scale is very tricky and should be clearly justified in each 
application. For instance, for traffic arrival curves, it could 
be set as the time until buffer overflows, or it could be 
set as a value that has been verified from real trace data. 
Nevertheless, as we have seen from the above example, 
the bounding functions are very sensitive to the maximum 
time scale. Never use a small maximum time scale value 
for the purpose of obtaining a tight performance bound 
without clear evidence demonstrating why the value is 
practical for the application in consideration. 

• Rule 3: After model transform, another constraint should 
be posed on the selection of 9, that is, the value of 9 must 
make the system stable. Specifically, after the transform 
the value of h(a, 0) must be bounded, where a and (3 
are the arrival curve and the service curve, respectively. 

Following the above rules, we can obtain tighter bounds 
shown in Figs. [3] and [4] if we directly build a v.b.c model using 
existing results from queueing theory O, iTPSl . specifically, the 
result on the (virtual) queueing length distribution with time 
t. 

C. Further Discussion: 9-m.b.c. Model And 9-s.c. Model 

In the above guidelines, the most difficult part is to handle 
the time-increasing problem. Another method is proposed to 



avoid the time increase problem in [ 1 1 1 by using a new traffic 
arrival model and a new service model, namely 9-m.b.c model 
and 9-s.c. model. We only discuss the 9-m.b.c. model, since 
the same conclusion is applicable to the 9-s.c. model. 

Definition 10: ifTTI A flow is said to have a 9-m.b.c. 
stochastic arrival curve a G J- with respect to 9, with bounding 
function f e (x) £ T, if for alH > and all x > 0, there holds 
for some 9 > 0, 

Prob{ sup [ sup A(fi, s)—a(s—fj)—d-(t—s)] > x} < / (x). 

0<s<t 0</j.<s 

(39) 

Clearly, the 9-m.b.c. model is a scaling-up method, i.e., it 
raises the curve a to get a tighter bounding function, hopefully 
to obtain a bounding function that is independent of time t. 
It has been proved lUTl that if a flow has a v.b.c stochastic 
arrival curve a with bounding function f(x) G Q, it has a 9- 
m.b.c stochastic arrival curve ag with bounding function f e , 
where for any 9 > 0, 



/(*) + \ I f(v)dy 



a e {t) = a{t) + e-t,f[x) 



• Rule 4: Generally speaking, if we can model a flow with 
the 9-m.b.c. traffic arrival curve such that the bounding 
function is independent of time t, we can find for the 
flow a v.b.c. traffic arrival model with a bounding func- 
tion independent of time t as well. Nevertheless, if the 
time-increasing problem exists in a weak model, model 
transform to a 9-m.b.c. model may not be very helpful. 

VIII. Conclusion 

An useful analytical technique for performance evaluation, 
stochastic network calculus has not got the deserved fame and 
its application in practice is lagging far behind the theoretical 
development. This should not be the case. This paper is to 
provide the guidance in the practical use of stochastic network 
calculus. Following the suggestions and the two-step approach 
in the paper should help the novice quickly grasp this useful 
technique. To conclude, we borrow Kleinrock's last words in 
his classical queueing theory book lfl4l : It now remains for 
you, the reader, to sharpen and apply the new set of tools. 
The world awaits and you must serve! 
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Appendix 
Proofs of Results 

Proof of Lemma [101 Given a Poisson arrival process A(t) with 
a rate A, based on Chernoff bounds 1 18 1, we have for Vf? > 0, 

P{A(t) -Xt>x} = P{e eA{t) > e 0{xt+x) } < f (xt+x) > ■ (40) 

It is easy to see that Ete ^) = e At ( e£ '- 1 ' (e.g., refer to Lemma 
5.3 in 1181 ). Therefore, «40t is equivalent to 



P{A(t) - Xt > x} < 



\t(e»-l) 



e B(\t+x) ' 

To obtain the tightest bound of the above inequality, we calculate 



(41) 



o At(e"-l) 
,6(\t+x) ■ 



= mm(At(e -l)-e(At+a;)) 



which can be obtained by calculating 



do 



From Equation <42L we obtain 



= In 



\t + x 
Xt ' 



(42) 



(43) 



Finally, the theorem is proved by applying ( 143b into d4U . 



